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LIST or SYMBOLS 


1 . Capital letters denote matrices unless otherwise stated. 

2 . Lov/er case letters denote column vectors unless other- 
wise stated or clear from context. 

3. A* denotes the matrix conjugate transpose of A . 

4 . A~* denotes matrix inverse for nonsingular A . 

5. A + denotes the generalized inverse of A . 

6. H will denote a hermitian idempotent matrix (h.i.) 

i.e. a matrix such that H* = H and HH ■* H . 

« 

7 . R(A) denotes the range space of A i.e., the collec- 
tion of all images of column vectors under the trans- 
formation A . 

8. P R /*\ will denote the orthogonal projection on the 
range of A i.e. a hermitian idempotent leaving R(A) 
fixed . 

9. E n will denote m - dimensional euclidean space. 

10. diagCa^ a 2 , ...» a n ) denotes a diagonal matrix. 


/ 


INTRODUCTION 


The primary concern of this paper is to Investigate 
the problem of Inversion of singular or non-square matrices 
In this connection, a new algorithm for computing the gen- 
eralized inverse of an arbitrary complex matrix is given. 
For a non-singular matrix the algorithm gives the ordinary 
inverse of the matrix. 

The paper is divided into several sections. The 
first two sections give a definition-theorem expose' of 
the known results in the literature. The following section 
give a new explicit form together with an algorithm for com 
putirig the new explicit form. An application to least 
squares approximation Is given that can easily be realized 
in trajectory analysis problems. Finally, a computer pro- 
gram for computing the generalized inverse of a matrix is 
given utilizing the algorithm mentioned in the latter para- 
graph. 
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DEFIMI T10HS AMD EQUIV A LENT F ORMS 

A. BJerhammer [2] 1 , E. H. Moore [10], and R . Penrose 
[}1] independently generalized the concept of matrix in- 
version to include arbitrary complex matrices. The gen- 
eralized inverse of a singular or non-square matrix possesses 
properties which make it a central concept in matrix theory. 

We will give a definition-theorem expose, inserting 
where applicable, references and special problems. The 
following fundamental theorem due to Penrose [11] wi]l be 
stated without proof. 


THEOREM 1. The four equations 


(1) AXA - A 

(2) XAX * X 

(3) (AX)* - AX 

(4) (XA)# * XA 


have a unique solution X for each complex matrix A . 


Definition 1. The solution X in. THEOREM 1. will be de- 
noted X = A + and called the generalized inverse of A . 


The following theorem gives an equivalent form of A 


THEOREM 2. For any mxm matrix A over the complex 
field, X « A + is the unique solution to the equations 

lumbers in brackets refer to correspondingly numbered 
papers in the references. 


O 
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AJt c P R(A) 

and 

XA “ P R(X) 

where R(A) is the range space of A in E m and 
is the orthogonal projection on R(A) . 

Froof: THEOREM 1. implies that AX is a hermit i an idern- 

potcrit (see list of symbols) leaving A fixed i.e., 

(AX)A c A . Hence AX must be a projection. We mav 
conclude the same about XA . 

V/e proceed to give properties of the generalized in- 
verse and possible computing schemes. 

THEOREM 3 . Let A be an arbitrary complex matrix. Then, 


for 

scalar 

X f 0 and unitary U and 

(a) 

A + (A + ) S 

A* = A + = A*(A + )*A + 

(b) 

A + AA li * 

A* = A*AA + 

(c) 

(A + ) + = 

A 

(d) 

(A* ) + = 

(A + ) * 

(e) 

A + = A' 

for nonsingular A . 

(f) 

(XA) + = 

>A + 



X 

(E) 

(A*A) + 

= A + (A + )» 

•(h) 

(UAV) + 

« V -1 A + U - ^ 


/ 
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(1) 

V D* 

and 

A i A j 


Vi " 0 

for 

i / J 


imply 

A 

*■ E 

A i 

(J) 

If A is 

normal (i 

. # 

.e . A 


then, A + A 

* AA + and (A n 

(k) 

A, A A, A + 

and A + A 

all 


to trace 

(A + A ) . 


(X) 

A = (A A) 

+ A* 



We note that ( 1 ) reduces the problem of computing A + 

to that of computing the generalized Inverse of a hermitian 

* 

matrix A A . Moreover, such a matrix can always be 
diagonalized by a unitary transformation i.e., 


D = U(A A ) V - diag (a 1 , . . . a n > 


Nov; (f) and (h) Imply 


(A*'a) + ■= VD + U = V diag (5 £ n )U 

V/e tacitly assume that If a^ = 0 then the corresponding 
term In diag (^. , 1 ) is zero. It is not usually an 

easy task to determine the unitary transformations U and 
V . Methods for computing the generalized inverse have 
geen given by various authors [2], [3], [7], [8], [12]. 

The following is a theorem of major importance charac- 
terizing all solutions of the matrix equations AXB = C 


which have some solution X . 


o 
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THEOREM . For the matrix equation AXB ** C to have a 
solution, a necessary and sufficient condition is: 

AA + CB + B « C 

in which case the general solution is: 

X ■= A CB + + Y - A + AYBB + , 

# 

where Y is arbitrary (to within the limits of being con- 
sistent with dimension in the indicated multiplications) [31]. 

Proof: Suppose X satisfies AXB e C . Then, 

C * AXB ■= AA + AXBB + B = AA + CB + B 

''inversely. If C » AA + CB + B then A + CB^ is a particular 
so.dtion. Clearly, for the general solution we must solve 
AXB = 0 . Any expression of the form 

X = Y - A + AYBB + 

is such a solution. Moreover, if AXB « 0 then, 

X = X - A + AXBB + 

0 

We note that the only property required of A + and B + in 
the theorem is AA + A *= A, BB + B * B . 

COROLLARY 1. The general solution to the vector equation 

Px = c is 

x «= P + C+ (I - P + P)y . 


O 


V 
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where y is arbitrary, provided a solution exists. 

COROLLARY 2. A necessary and sufficient condition for the 
equations 

AX * C 

and- 

XB ' D 


to have a common solution is that each have a solution and 

AD *= CB 


Proof: If AX * C and XB - D have a common solution 

then clearly each has a solution and 


so that 


AXB = CB 
AXB = AD 

CB * AD 


> 


In order to obtain the sufficiency we set 

X *= A + C + DB + - A + ADB + 

which is a solution if AD * CB, AA + C * C, and D3 + B = D . 


THEOREM 


5. We have: 

(1) A + A, AA + , I-A + A, and I - AA + are h.i. 
(see list of symbols) 

(2) H Is h.i. inplios H + » H 
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Proof: The proof requires a straightforward application 

of THEOREM 1. 

In general, the reversal rule (i.e., (AB) 4 * B 4 A + as 
in the case of the standard inverse) does not hold. R. Cline 
[5] recently obtained the following result. 

THEOREM 6. Let A and B be matrices with the product AB 
defined. Then, 

(AB) + « 

where AB * Aj^ 

and Bj * A + AB 

A 1 c AB i B i 
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THK EXPLICIT FOR" 

• 

utilizing the properties of A + in the proceeding 
sections, we develop an explicit form which gives rise to 
an algorithm for computing the generalized inverse of an 
arbitrary complex matrix [7]. 

THEOREM 8. For any matrix A, A + *= WAY, where W and Y 
are any solutions of 

(1) WAA *= A* 

and 

(2) A AY « A 


Proof: Equations (1) and (2) indeed have a solution 

V/ *= y r= a + . Moreover, if W and Y arc any solutions v:e 

have 

AWAA* ” AA and A*AYA = A A 

so that AWA = A and AYA *= A 

(Note: BAA* = CAA* implies BA = CA). 

In addition, 

x s x * * x r. * 

WAA W = A W and Y A AY = Y A 

imply (V/A ) * = WA and (AY)* *= AY 


If we let X = WAY, X satisfies the four equations of 


so that A = 




THEOREM 1 


X = WAY. 
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4 * * 

COROLLARY 2. For any matrix A, A *= A v/here 

and are, respectively, any solutions of 


(AA )S (AA ) c (AA ) 

and 

(A*A)S 2 (A*A) - (A* A ) 

* 

Proof: According to THEOREM 3. we have that W = A Sj 

£ 

and Y *S ? A are solutions of equations (1) and (2) of 
THEOREM 3. provided 


(AA*)S 1 (AA*) « (AA*) 

and 

(A*A)S ? (A*A) = (A A) 


The corollary follov/s. 

THEOREM 9. If B is a matrix and there exist nonsingular 
matrices P and Q such that PBQ = E is an idempotcnt 
then B « QEP is a solution of BXB = B . 


Proof: 

If P, Q, and 
_1 

E satisfy the hypothesis 

-1 

of the 

theorem 

then B = P EQ 

and 



B3B = (P _1 

EQ“ 1 )QEP(P~ 1 EC! -1 ) = P" ^EQ 

~ 1 * B 


COROLLARY 2. and THEOREM 9. suggest an algorithm for 
computing the generalized inverse of a complex matrix F. 

4 - £ 4 - * 

Consider the equation F « (F F) F [6], which reduces 

the problem of finding F + to that cl finding the general- 

* 

lzed inverse of the hermitian matrix F F = C. Since 


:o 


p * 2 

(C ) * C , there exist nonsingular matrices P and Q 

(products of elementary matrices obtained by simple 
elimination) such that 



I 


o 


where I is a rank r identity matrix and the Z are 
r ' 

zero matrices. We set C *= A In COROLLARY so that 

* * * n 2 

A A r AA =CC=CC'=C. According to THEOREM 9« 
choose solutions e QI q P so that C ~ (CS^) C, 

(F*F) + *= C + , and finally, 


Computing programs for calculating S^ and S^ are 
now in existence (e.g., STORM, Statistically Oriented 
Matrix Program, IBM). In general, these programs only 
compute some solution of the equation AXA = A, usually 
different from A + . These results allow one to construct 
a solution to all four Penrose equations (THEOREM 1.), 
given only a solution of the first, namely, AXA = A . 


APPLICATION TO LEAST SQUARES APPROXIMATION 


We will now state an application that can be realized 
in trajectory analysis problems. For the sake of simplicity 
we will not consider weighting and only mention that weight- 
ing introduces no difficulty. 

The vector equation Ax *= b does not, in general, 

4 

have a solution x . However, all candidates for a least 

squares solution (i.e., a solution vector x minimizing 

* 

(Ax - b) (Ax - b) must be solutions of the normal equations 

A*Ax *= A*b [8] 

THEOREM 7. Let A be any matrix (mxn) and b be any 
vector (mxl) . The equation 

* # 

A Ax ■ A b 

always has a solution and hence a general solution given by: 

X *> (A*A) + A*b 4- (I - (A*A) + A*A)y 
= A + b + (I - A + A)y 

x 

Moreover, if A A is non-singular then the solution is 

x *= A + b 

and is unique. 


Proof: 


We will first show that 



1 ? 


has a solution. Consider the vector: 

x * A + b 

Lince THEOREM 3 • ( b ) implies A*A(A + b) *= A*b vre have that 
x * A + b is indeed a solution of (1). The existence of 
this solution together with COROLLARY 1. implies that the 
general solution to (1) is: 

(2) x « (A*A)V'b + (I - (A*A) + A*A)y 

Using THEOREM 3.(1) we see that 

x *= A + b + (I - A + A)y. 

Finally if A A is non-singular then 

x ■= (A*A) + A + t + (I - I )y 
«= A + b 

and (1) has a unique solution. 

In summary, we know that if x is a least squares solution 

of Ax = b, then x must satisfv 

* * 

A Ax 5 A b 

All solutions of this equation are given by x * A + b + (I - A**A)y. 
Any vector of the form 

x “ A + b + (X - A + A)y 

is a "candidate" for a least squares solution and this form 
describes the "class of all candidates.” 
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COROLLARY 3. Every solution of A Ax * A b minimizes 

f 

q e (Ax - b) (Ax - b) provided Q has a minimum. 

Proof: V/e know that any vector at which Q is minimum is 
of the form 

x ■ A b + (I - A + A)y 

If Q has a minimum let 

x-j * A + b + (I - A + A) y ? 

be any other solution. We will show that 

(Ax^ - b)*(Ax^ - b) ** ( AXp - b)*(Ax 2 - b) 

To do this we examine Ax^ and A:< 2 using THEOREM 1. 

Ax-^ ~ A(A + b + (I - A + A)y^) 

*= AA + b + (A - AA + A)y . 

■= AA + b + (I-Dy, 

■= AA + b 

Similarly Ax_ = AA + b 
so that 

(Ax j - b)*(Ax 1 - b) = (Ax 2 - b)*(Ax 2 - b) 
that is, every vector of the form 

x = A + b + (I - A 4 A )y 


yields the same minimum value of Q . 
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SUBROUTINE OEN1 NV 

GENINV is a FORTRAN IV subroutine, written by L.F. 
Guscnan, Jr., Theory and Analysis Office, which is used to 
compute the generalized inverse of an mxn matrix A . 

t 

All computations are done in double-precision floating point 
arithmetic. The subroutine is based on the algorithm sug- 
gested by the explicit form. 

CALLING S EQUENCE 

CALL GENINV (A, AP, M, M, L. E), 

where, 

A is a double-dimensioned , double-precision array 

containing the original matrix. A is dimensioned 
A (25 , 25). 

AP is a double-dimensioned, double-precision array 

where the generalized inverse of A will be 
computed. AP is dimensioned AP(25, 25). 

M is the number of rows in the original matrix. 


N is 

the number of columns 

in the original matrix. 

L is 

tv; ice N. 


E is 

some small number for 

near-zero divisor test. 

METHOD 



Given A 

• 

(PRINT A) 

Compute: 

C = A* A 

(PRINT C) 


C 2 *= cc 

(PRINT C ? ) 


O - 15 - 

Find non-singular matrices F. and P such that 



PRINT E, P, I Q 

(A form of Gaussian elimination with pivoting employed) 


Compute: 

R *= PI q E (PRINT R) 

then 

C + ■= CFtCRC (PRINT C + ) 

also 

A + ■= C + A* (PRINT A + ) 

o 

Remarks 

The program uses two double-precision arrays 
CSQ(50, 50) and B(25, 25) for internal manipulation. The 
subroutine leaves the original matrix A intact. 

Results are printed after each step as indicated. 




O 
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SU3KCUI I N £ 0 L N l N V ( A , A P , M , N , L , C ) 

DIMENSION A ( 25,25) t API 25, 251 , C SO { 50, 50 ) , H ( 2 5 , 25 ) 


c» 

DOUBLE PRECIS 1UN C SO , R , l J I GA , X , D A C S , A , A P 


c«*« 
c • • • * 

THIS SLJROUTINE COMPUTES THE GENERAL I ZtO INVERSE 
algorithm fiy h. r. dlceh , 

OF A MATRIX 

c« • 
c«* 

CALLING SECUtNCE 


c * * 

c*« 

CALL GtNlNV( A,Ap,M,N,L,E ) 


c* * 

c» * 

A ( M , N ) - LOCATION OF ORIGINAL MATRIX 
A P ( N , M ) - LOCATION Of COMPUTED GENERALIZED INVERSE 

C« * 

c«« 

K - NUMBER Of ROW S IN ORIGINAL MATRIX 
N - NUMBER OF COLUMNS IN OK If, INAL MATRIX 


c «• 
c* * 

L - F*N 

E - SMALL POSITIVE NUMBER EUR NEAR-ZERO DIVISOR 

TEST 

c* 

c 

INITIALIZATION 


N P 1 = i\ + 1 
KR l T E ( 6 , 1 00 ) 

100 

FORMAT ( 1 H 1 , l 3H THE MATRIX A//) 
DO 5C I = i , y 


50 

WRITE! 6,200) ( A € I , J) , J = 1 , N) 

continue 


c 

c 

* 

COMPUTATION OF C = A A 



DO l 1=1, N 
1 1 - I + N 


CO 1 J * 1 » N 
J 1 = J ■» N 

C SQ ( II, JD-O.OUO 
DO 1 K = l , M 

1 

CS(J( I 1 , J1 ) =CS'J( I 1 , J 1 ) 4- A (K , I ) *A(K, J ) 
CONT INtE 


10 i 

W R I T E ( 6 , 101) 

FORMAT ( 1H0, l 3H THE MATRIX C / / ) 


DO 5 1 I =NP 1 , L 

KRI TE( 6,200 (CSQ( I , J) , J = NPi , L ) 

51 

200 

CUN I I- ILL 

FORMAT (1H ,6021.12) 


C 

c 

2 

COMPUTATION OF C = CC 

0 

DO 2 1=1, N 
11 = 1 +.M 


DO 2 J= 1 ,N 
J1=J4>J 

• 

CSO( I » J ) =0.0C0 
DO 2 K = 1 , N 




no! I ! o o o o o; 


O 
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K I *= K ♦ >4 ~ 

CSQ< 1 1 J)*C$C( I t J) ♦ CSUl I IjKl ) t JD 

2 CONI I 4UC 

V.'R ! T t ( 6 » 1 0 / ) 

107 FORMAT tlH0 9 I ?H THE *•'. A T R I X CSQUARL//) 

D0_5 IJ) I = 1 » M 

NR l TEC 6»200) (CSOU.J).J-lt N ) 

510 CUM I ML 

C 

C DUILD augment ed matrices 

DO 5 1.1,N 

DO 5 J=\P1 , 1 

IF C ( j-.\ J - I ) 3 ,4,3 

3 CSQI It J)«O. OCO 
GO T l 9 

A CSJ( 1, J)«U OQO 

5 CONTINUE 

_D0_8 J = 1 i N 

DU 8 I = .\P 1 , L " 

I F { ( J)6, 7 f 4 

6 CS'M I , J ) =0 • Oo 0 
CO To P 

7 CSQt If J) = 1.000' 

8 CONTl'MUt 


C QMPtT AT ION Of I 

0 

A FORM OF GAUSS I AN I L I M I NA I I ON IS EMPLOYED 

_ DC) 2 2 K = I , N 

KP 1 = K ♦ 1 

IR = K 

JC = K 

B1GA«0AQ S(C SQ(K,K) ) 

DO 10 r.KtN 

DO 1C J=K,N 

I F ( 9 I G A - C 53 S ( C $ N ( It J)M9 9 10»10 
9 IR«I 
JC - J 

BIGA.OADSICSQ ( I , J ) ) 

10 CONf I nUE 

I F ( B IGA-E ) 2 3 f 23,11 


EXCHANGE ROWS 

11 IK fR-K 112,14, 12 

12 DO 13 J - 1 * L 
X=CSL ( IK , J ) 

C SQJ IR, J)=CSC(K,J ) 

CSQ ( K , J ) >X 

13 CONT I IUE 
C 

C EXCHANGE COLUMNS 

14 IF'(JC-K)15, 17,15 

15 CO 16 I = l , L 
X=CSG( l , JC ) 

CSQJJt JC)=CSG( I,K) 


0 


13 — 






CSC I I • K » «X 

16 CONTINUE 

c 

C DIVICF ROW K BY CSC(K,K) 

17 X=CSu(K|K) 

DO 18 J*K,l 

CSOUt J r«cs"6 (K» J)/x“ 

18 CONI I ILL 

k ik-nI i WiTZti'2 

c 

c zero cclumn k crio* the diagonal 

19 00 ?C I «KP 1 t N 

XrCS'. ( I t K | 

DO 2C J = K » l 

CsoTit jT^c sq n ,"jl-x»cio JkTj i 

20 CONTIN UE 

c 

C ZERO ROW K TO THE RIGHT UF THE 01 AGONAL 

DO 21 J - K P 1 » N 
X*CSCI K|J) 

00 2 1 l - L » L 

CSOI I » J)*CSQ( I » J ) - A * C SQ ( I fK| 

2 1 CONTI iUE 

22 CONTINUE 

2 3 C U N T l 1 Ut 
WRITE! 6 > I02j 

102 F0RMAfriHC7l7H T'h! H a T K I X I/ l m.//) 

DO 52 1*1 tN 

WR I T L I 6 , 200 ) ( C S Q ( I » J ) « J = 1 » \ ) 

5 2 c orn uu 

WRlTt (6, 105) 

105 FORMAT ( lHO > i 3H THE MATRIX f / / ) 

CO 5C2 IVT.N 

WRITE(6t200) ( CSQII t J) , J«NP1 tL) 

502 CO O I |UE 

ViR I T £ ( 6, 106) 

106 FORMAT ( lhO * 1 TH THE MATRIX P//) 

DC 5 C 3 1»N° 1 1 L 

WRITE(6#2G(T) ( C SC ( l ,J) , J « 1 , N ) 

503 COM MU 

C 

c computation of r=pi e 

c “ o 

00 2 A I »ltN 
DO 2 A J*l f N 

B( I , J) =0.000 

DU 2A K*l«N 

J 1 = J*N ; 

B ( I » J ) "5 ( I t J ) +CSQ ( I » n ) *CSQ ( K , J 1 ) 

24 CONTINUE 

00 25 1*1 f N 

Cu 2 5 J- 1 |N 

C SQ ( I» J) =0.000" 

DO 25 K = 1 * N 


G 


11 = 1 4N 


► 19 - 

25 COST l.MLt 
WRIHI6, ICO) 

103 FOftMAtl iHltlTH T hi MATRIX k//) 

CO 5J0 I ' 1 ♦ 'J 

WRlT€(6t200MCSQ(l t J) , J- l ,.\) 

_ 530 cimnuL 

c ♦ 

C COMP UTATION Uf- C •£R(R£ 

CO 2c |«1 
DO 2t J = 1 1 N 
Jl= J4N 

e< i , j)rQ.oco 

00 26 mTUT 

Kl«K 4,N 

B ( I » J ) «0 ( I V J ) + C SQ ( l • K ) *CSC ( K l , J 1 ) 

26 COM I viuf 
CO 27 i»l f N 

CO 2 7 J g 1 >N 

CSQ( I, J)«0.000 

op 27 K«I»N 

CS'j( I , J) vCS'M I , J ) 451 I , ) • ( K, J ) 

27 CON I I iUL 

DO 2d I*l#N 

11«I iH 

OU 2 3 J - 1 , N 
B( I , J) -0.000 
DO ?t K«1»M 
K 1 = K ♦ . I 

B( I , J ) I , J) 4CS 3 ( I 1 «K1 ) •cscfl K, j ) 

2fl C ONT I NUE 

WRITE (6 » 1031 

103 F ORMAT C IH Qt 1 7H THE MATRIX CPLUS//) 
DO 53 I « 1 , N 
WRIT l ( 5,200) (Btl f J) t 
53 CONT I IOC 

C ♦ 

c ' computation of a =c a 

OU 25 |*l tM 
00 ? J = I » F 
AP { I , J ) -C. 003 
00 25 K«lfN 

A Pi I , J ) = A P ( It J)+D( I ,K) «_A ( J,K) _ 

29 CONTINUE 

W RIT£(6t 10 4) 

105 FORMAT ( IhO , 17H THE MATRIX APLUS//) 
CO 5A I=l,i 

I T L ( 6 .200 ) ( Af’ (I , J) 

_ 5 4 CONT I ^UE 

RETURN 

END 



o 


o 


n 


THE MATRIX A 


0.4COOOCOCCCOOD 

01 

-0. 1CCCOCO0OOOOD 

01 

-0.3000000000000 

01 

0.2000009000003 

01 

-0.2C000CCCC00UU 

0 1 

0.500000O00000D 

01 

-o. 1000000000000 

01 

-c .^oooocooooooo 

01 

0.2CCCCCCC0GGCD 

01 

0. 30000000CCCCD 


-0.9000000000000 

01 

-C. 5000000000000 

01 

the matrix c 








0.24Q0CCCCCCCQD 

0? 

-o. eoccocoocoooo 

01 

-0.280000000000D 

02 

0 .4000000000000 

01 

-0.8C0CCCCCC000D 

01 

0. 35000000CCCOJ 

0 2 

-0.29C000000C00D 

02 

-0. 3200000000000 

02 

-0.28000C00OCCCD 

02 

-0.29000C0CC000D 

02 

0.9100000000000 

02 

0. 42 000 00 000 CO J 

02 

0. 4CC0CC0CC00C0 

01 

-0. 320000C000000 

0? 

0.4200000000000 

02 

0.380000000000D 

02 

THE MATRIX CSCUAR' 

c 







0. 144CCC0CCCC0D 

04 

0.212000000000D 

03 

-0.2020000000000 

04 

-0.6720000000000 

03 

0. 21200 OOCOOCOO 

03 

0.3154C0CGC000U 

04 

-0.47 74000000000 

04 

-0.3586000000000 

04 

-0. 2H200C0CC00GD 

04 

-0.4774OC000C00O 

04 

0. 1 1 67000000000 

05 

0.62340000C00CD 

04 

-0.6 /2CGCCCCC000 

03 

-0.3086000000000 

04 

0. 6234000000000 

04 

0.4248000000000 

04 

THE MATRIX I ZERO 








l.OOOCOCOCOOOOD 

00 

-0.2776557561560- 

■16 

0.555111512313D- 

• 

16 

-0.1387778780780- 

•16 

0. 


l.OOOOOOOCOOCOO 

00 

-0.5551115123130- 

16 

-0. *-551115123130- 

16 

-0. 


c. 


l.OOOOOOOOOCOOD 

00 

0.55511 15123130- 

16 

0. 


0. 


-0. 


0. 4209965 709 3 CO- 

12 

the matrix e 

0 . 


0. 


0.8568900291350- 

C4 

0 . 


0. 


0 . 8326 I 36 8 4944 D- 

■0 3 

0.34 06 08 20 32 5 CO- 

03 

0. 


0. 


0.3505 146 32 5360- 

•01 

— 0.7372689575180- 

02 

0.4064405973610- 

01 

1 .OOOOCCCCOOOOD 

00 

0.3774753283730- 

14 

0. 7272727272730 

00 

-0.9090909090910 

00 

THE MATRIX P 

0 . 


0. 


0 . 


1.0000000C00009 

00 

0. 


1.0000000000000 

00 

0.8624011351510 

00 

-0.6494804694060- 

14 

1 .ocooccccooooo 

00 

0.4090631 191 090- 

•00 

-0. 1813964850720- 

00 

0.7272727272739 

CO 

0. 


u. 


1 .0000000000000 

00 

-0.909090 V09U9 10 

CO 


o 


THE MATRIX A* APLUS 

l.COOOOGCCOCOOU CO 

-0.5953570969550-14 

-0.2745026428390-13 

• 


0. 5d2867087928D- L4 
-0, 1265t>54?4eC7Q-l 3 

1 .0000000000000 00 
0.9159339953160-14 

0. 69388939039 10-14 
1.0000000000000 00 



THE MATH T X A*APLUS*A 





0.4C00CCCCCC00D O’ 

-1 .COOCOCCOOOOOD 00 

-0.3000000000000 01 

0 • 200000000000D 01 


-0.2C000C0CCCG0IJ 01 
0.200C0C0C0G000 01 

0.5000000000000 01 
0.30000000CCC00 01 

-l.OOOOOOOOOOOCD 00 
-0.9000000000000 01 

-0.3000000000000 01 
-0.5000000000000 01 


THE MATRIX APLUS* A*APLUS - 




0. 1922H07C I L7540-00 

0.6035087 719300-01 

-0. 3649122807020-01 

• 

. 


0.2COCOCOOOOOOD-CO 
-0 . 56 1 AO 350 8 7 72 0-02 

0. 3 OOOOCOOOCCOU-OO 
0. 5298245614040-01 

-0.1000000000000-00 
— 0.90175438 59 6 50— 01 



0.20 70 i 754 3860U-CO 

0. 1087 719298250-00 

-0. 1122807017540-00 



THE MATRIX A P L 'J S * A 

0.5754385964910 CO 
0.6439293042830- 14 

-0.2534004053710-13 

l .ooooooccoocoo on 

-C« 30877 19 298 2 50-00 
0.3202993426040- 1 3 

0. 3359649122810-00 
0. 328070903777D-13 

• 

-0. 3087719293250-C0 
0.38596 4 * 9122 tt In- 00 

-0.7521760991 «4u-l4 
—0.286992651 8660—13 

0.7754385964910 00 
0.280 7017543860-00 

0.2807017543860-00 
0.64912280701 80 CO 



THE MATRIX R 

0.42099657093^0- 12 0, i 589 1 6G29357j> 

-0.76382 71 398590-17 0.3 106 105263160 

0*3061 3 2668 3280^12 -0*601760 7655500 

-0.3826686.44247D-12~ 0* 35051 4832 S3~6f> 

THE MATRIX CPlUS 

0 . A 194 5 706 37120-01 0 . 6 02 1 0 52 63 1 5 80* 

0.602 1052631580-01 0.1400000000000- 

0.54 086 79 593 720-02 ’ 0.2 3 73947 363420* 

0* 50467220663 30-01 0*83263157 09470 

THE MATRIX APLUS 

0. 19228070 l 7540- CO 0. 603 50 8771 9300 

0.2C0OOCCCCCCCO-0O 0.300COCOOOOOQD 
■0. 5614035 087720-02 0.5 298 2 4 5614040* 

0 . 2 0 7 C 1. 7 5 4 3 6 6 0 0- CO 0. 1067719298250* 


Q 


26 0.306179^243130-12 -0 . 38 2 7? 4 l 55 3980- 1 2 
01 -0.601760765^500-02 0.3505148325360-01 
0? __ 0* 1562406843780—0 2 -0 *7 37268957546 0-02 
01" -0. 7372689575460-02 0 * 406 44 0 5 97365 D- u I 


01 0*5408679593720-02 0*5046722068330-01 

0 0 0*2378947368420-01 JO *852631 5769470-01 

01 0.1097026777470-01 0. 147257o 177290-01 

01 0. 14725761 77290-01 0 . 672 94 5 5 2 1 699[)-0 1 


01 -0.3649122 80 7020-01 

00 -0. lOOOCOOOOOOOD-OO 

01 -0.9017543359650-01 

CO -0*1122307017540-00 
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